EXACT SOLUTIONS FOR WAVE PROPAGATION IN 
BIREFRINGENT OPTICAL FIBERS 



> 



o 



E. Alfinito, M. Leo, R.A. Leo, G. Soliani and L. Solombrino 

Dipartimento di Fisica dell'Universitd, 73100 Lecce, Italy 
and Istituto Nazionale di Fisica Nucleare, Sezione di Lecce. 



O 

, Abstract 
(N ■ 

We carry out a group-theoretical study of the pair of nonhnear Schrodinger 
equations describing the propagation of waves in nonhnear birefringent optical 
fibers. We exploit the symmetry algebra associated with these equations to 
' provide examples of specific exact solutions. Among them, we obtain the soliton 

■ profile, which is related to the coordinate translations and the constant change 

of phase. 



1 Introduction 

Q , The propagation of optical pulses in nonlinear birefringent fibers is described by the 

pair of nonlinear Schrodinger equations ^ 

Ai = i Ux + Utt + kv + {a\u\'^ + P\v\'^)u = 0, (1.1a) 

A2 = ivx + vtt + ku + {a\v\^ + p\u\'^)v = 0, (1.16) 

where u = u{x,t) and v = v{x,t) are the circularly polarized components of the 
optical field, x and t denote the (normalized) longitudinal coordinate of the fiber and 
the time variable, respectively, k is the birefringence parameter and the coefficients a 
and f3 are responsible for the nonlinear properties of the fiber 0, |^, ^. Performing 
the change of variables x — > ax, t — > /it, with jj? = ±|a, Eqs.(l.l) take the 
form 

iUx ± -Utt + k'v + (1^1^ + (T|f 1^)^ = 0, (1.2a) 



i ± -vtt + k'u + {\v\^ + a\u\^)v = 0, (1.26) 

where k' = k/a, + (-) holds in the anomalous (normal) dispersion regime and c" = f = 
jz^, B being the third-order susceptibility coefficient H H- For a = (3 (a = 1) 
and = 0, the system (1.1) has an infinite set of constants of motion and may be 
solved by the inverse scattering method 0. 

In general, i.e. for a ^ (3 and 7^ 0, the system (1.1) is not integrable by 
inverse scattering. In this case, Eqs.(l.l) possess three constants of motion only 
[|l|. On the other hand, it is well-known that a powerful tool for handling both 
integrable and nonintegrable differential equations is represented by the so-called 
symmetry approach . This method, which is based on the Lie group theory, consists 
essentially in looking for symmetry transformations that reduce the equations under 
consideration to certain ordinary differential equations, each of them comes from an 
invariant quantity associated with a given symmetry allowed by the system. Following 
this idea, in this work we apply the symmetry approach to Eqs.(l.l). We display 
examples of new exact solutions in both the cases a = (3 and a ^ (3. In this regard, 
we observe that the birefringent parameter k involved in Eqs.(l.l) is real. However, 
the symmetry algebra found for a = /? does not depend on k. This fact has suggested 
us to study the system also for imaginary values of k. In such a situation, at least when 
u = V, k can be interpreted as the loss coefficient of the fiber [|^. Anyway, we have 
obtained an interesting exact solution for -u 7^ f as well. This solution is derived from 
the Galilean boost. Another important aspect of the symmetry reduction technique 
is the determination of the infinitesimal operator which is responsible for the soliton 
profile and the periodic configuration. 

In Sec. 2 we outline the method of symmetry reduction and obtain the symmetry 
algebra and the corresponding symmetry group related to Eqs.(l.l). Sec. 3 contains 
examples of specific exact solutions, while in Sec. 4 some concluding remarks are 
reported. 

2 The method of symmetry reduction 



2 



2.1 a) The symmetry algebra 

The method of symmetry reduction (SR) consists of an apphcation of the Lie group 
theory to reduce Eqs. (1-1) to a system of ordinary differential equations. 

A fundamental step of the SR procedure is to obtain the Lie point symmetries 
of Eqs. (1.1): in other words, the symmetry algebra L and the corresponding 
symmetry group Q of the equations under investigation. Then, we can build up 
solutions that are invariants under some specific subgroup of Q. The SR can be 
carried out via the determination of the invariants of Qq. Invariants are furnished by 
the partial differential equations 

V,/(a;,t,M,t;,n*,0 =0, j = l,2,...,n, (2.1) 

where {l^ } is a basis of the Lie algebra £o of and n is the number of the independent 
elements ( infinitesimal operators) Vj of £o- Once the invariants related to are 
known, Eqs. (1.1) can be written in terms of them. In such a way, we are led to a 
set of reduced equations which may yield exact solutions to the original system (1.1). 
The Lie point symmetries of Eqs. (1.1) can be found by resorting to the standard 
technique outlined in 0. Precisely, let us introduce the vector field 

y = G^x + ^2dt + ^3du + iidu^ + ^5^. + Ud,. , (2.2) 

where (j = 1,2,..., 6) are functions which depend in general on x, t, m, m*, f , f *, and 
'^x = -§^, and so on. A local group of transformations G is a symmetry group for 
Eqs. (1.1) if and only if 

j9r(2V[Aj] = 0, pA^^V[A*] = 0, j = l,2, (2.3) 

whenever = 0, A* = for every generator of G , where pr^'^'^V is the second 
prolongation of F . 

The conditions (2.3) constitute a set of constraints in the form of partial dif- 
ferential equations which enable us to obtain the coefficients The calculations 
have been performed in part by using the symbolic language REDUCE p|. We have 
achieved the following results. 
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Case I : a 7^ /3 

The symmetry algebra is defined by four elements, namely: 

Vi = d^, V2 = dt, V3 = i{udu - u*du* + vd^ - v*d^*), 

= xdt + -t [udu - u*du* + vdy - v*dy*), (2.4) 
where the nonvanishing commutation relations are 

[Vi,V,] = V2, [V2,V,] = ^Vs. (2.5) 

Case 11 : a = f3 

The symmetry algebra is of the sl(3,C) type. It is defined by eight elements : four 
of them coincide with the previous ones, while the others are given by 

Vs = xdx + ^tdt — ^{udu+u*du*+vdy+v*dy*)+ikx{vdu — v*du*+udy — u*dy*), (2.6a) 

= i{vdu — v*du* + udy — u*dy-), (2.66) 

Vj = {vdu + v*du* —udy —u*dy*) cos 2kx — i{udu ~ u* du* —vdy + v*dy*) sin2A;x, (2.6c) 

Vs = {vdu + v*du*—udy — u*dy*)sm2kx + i{udu — u*du*—vdy+v*dy*)cos2kx. {2.6d) 

The nonvanishing commutation relations fulfilled by Vi, . . . ,Vs, are (2.5) together 
with 

[V^i, V,] = Vi + kVe, [Vi, Vr] = -2kVs, [Vi, Vs] = 2kVr, 

[V2,V,] = ^V2, [V,,V4 = ^V,, (2.7) 

and 

[Ve, Vr] = 2Vs, [Vj, Vs] = 2V,, [Vg, V^] = 2V-j (2.8) 

The vector fields Vi,...,V8 are the generators of the infinitesimal symmetries trans- 
formations of Eqs. (1.1). These are the coordinate translations and the Galilean boost 
(V4), which are common to both the cases I and II. Moreover, for a = /? Eqs. (1.1) ad- 
mit the additional symmetry su(2,C) which is expressed by the generators Vg, V7, Vg 
satisfying the commutation rules (2.8)[|. 

^ It is noteworthy that for a ^ /3, the symmetry algebra turns out to be independent from the 
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2.2 b) The group transformations 

By integrating the infinitesimal operators Vi, Vg, we provide the group transfor- 
mations that leave Eqs.(l.l) invariant. These are, respectively: 

Vi : i — t, X — X + X, u — u, V — V, (2.9a) 

V2 : t — t + X, X — X, u — u, V — V, {2.9b) 

V^: x = x, i=t, u = u e*^, v = v e'^, (2.9c) 
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X 



X, i = t + Xx, u 



(At-iA^x) 



(2.9d) 




cos k{x — x) sin k{x — x) 
— smk{x — x) cosk{x — x) 




(2.9e) 



Ve : X — X, i — t, 




(2.9/) 
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t = t, 

_ giA(cos2fea; (T2— sin2fex as) 



U 



(2.9^) 



: X = X, t — t, 
— e 



^ \ _ „a(sin2fci; (T2+cos2fea; 0-3) ( 1 ' (2.9/i) 




parameter /s. Conversely, for a ~ j3 and A; such that Imfc 7^ 0, the symmetry algebra is defined 
by the vector fields Vi, V2, V3, and VJj. In other words, the presence in Eqs.(l.l) of a non-real 
coefficient in front of u and v changes the symmetry algebra, which is no longer of the sl(3,C) type. 
This fact is connected with the loss of the integrability property of Eqs.(l.l) for a = /3 and k such 
that Imfc ^ 0. 
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where A and di, (T2, are the group parameter and the Pauh matrices, respec- 
tively. In deriving (2.9), we have used the initial conditions x{X) |a=o= \x=o 
= t, u{X) |a=o = u, v{X) |a=o = V. 

Equations (2.9a)-(2.9h) tell us that iiu — f{x, t), v — g{x, t) is a solution of the 
system (1.1), so are 

^(1) = /(x- A,i), v^^^ ^ g{x-\t), (2.10a) 

= _ A), ^(2) =^(a;,t- A), (2.106) 
^(3) ^ e'^f{x, t), v^^^ = e'^g(x, t), (2.10c) 



u 



(4) 



/(x,t- Ax)et(^*-5^'^), 



,(4) 



g'(x, i — Ax)e2^ 2 



(At-iA^x) 



(2.10ci) 



u 



(5) 



= {/(e"^a;, e"^^) cos [A;x(e^ - 1)] + ig{e~^x, e~^t) sin [kx{e^ - 1)]}, 



^(5) ^ g-||^j.g-A^^ [/ca;(e^ - 1)] + if{e'^x, e'h) sin [kx{e^ - 1)]}, 



(2.10e) 



u 



(6) 



/(x, t) COS A + ig{x, t) sin A, 



= i/(x, i) sin A + g{x, t) cos A, 



(2.10/) 



' = (cos A — i sin A sin 2A;a;)/(a;, t) + (sin A cos2kx)g{x,t), 
v^"^^ — (—sin A cos2A;x)/(x, t) + (cosA + isinA sm2kx)g(x,t), (2.10^^) 



u^^^ = (cos A + i sin A cos2kx)f{x,t) + (sin A sin 2A;a;)5f(a;, t), 
v^^^ — (—sin A sin2A;x)/(x, i) + (cos A — isin A cos2kx)g{x,t). (2.10/i) 



3 Exact solutions 



As we have already mentioned, the method of symmetry reduction of a partial dif- 
ferential equation amounts essentially to finding the invariants ( symmetry variables) 
of a given subgroup of the symmetry group admitted by the equation under consid- 
eration. A basis set of invariants for the generators Vj can be obtained by solving 
Eq.(2.1). Alternatively, one can resort to a direct equivalent procedure by using the 
group transformations (2.9). 

The invariants can be exploited to provide exact solutions to Eqs.(l.l). By way 
of example, in this Section we shall deal with the invariants related to the symmetry 
operators a) Vq = Vi + V2 + V3, h) V4, c) Vi + V4 and d) V5. 

Case a) A set of invariants related to Vq is 

y — i — X — t — X, 

Ui = ue-'^ = ue-'\ Wi = ve''^ = ve-'\ (3.1) 
Inserting, for instance, 

u{x, t) = ^7(y)e"-, v{x, t) = W{y)e'\ (3.2) 

into Eqs.(l.l) (for a — (3) gives the pair of (ordinary) reduced equations 

iU' + U- U" -kW - a{\U\^ + \W\'^)U = 0, (3.3a) 

iW + W - W" -kU - q;(|C/|2 + \W\'^)W = 0, (3.36) 

where 

U^U{y), W = W{y), U'^^, and W^' = 

dy dy 

Now, let us look for solutions to Eqs.(3.3) of the type 

U{y)^p{y)e''^y, W ^ q{y)e''y (3.4) 
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where p, q are real functions of y and 7, 6 are real constants. In doing so, Eqs.(3.3) 
yield 

(27 — l)p' + kq sin {S — ^)y — 0, 



(25 - - kp sin (5 - 7)^/ = 0, 

p" + (7 — 7^ — l)p + kq cos (5 — 7)1/ + aip^ + = 0, 

q" + (5 - 5^ - 1)5 + A;pcos {5 - 7)y + a{p^ + = 0. 
, , dp , dq 



(3.5a) 



(3.56) 



Equations (3.5) produce some interesting configurations of Eqs.(l.l), such as the 
sohton profile and solutions expressed in terms of the elliptic Jacobi functions sn(-). 

To this aim, let us choose 7 = 5 = |. Then, by defining z = p + iq = pe''^^^ [p ~ 
\z\), we have 

'^'"^ ' " (3.6) 



where p' = ^ and c is a constant of integration. The sohton profile comes from (3.6) 



for c = 0, by taking a > 0, k < j. Precisely 



4 'a 



k{y- yo) 



(3.7) 



where yo is an arbitrary constant. With the help of (3.7), from (3.2) we obtain 

^i = -t; = et(*+^)p(t-x), (3.8) 

with pit — x) given by (3.7). 

We notice that for c = 0, and a < 0, k > j, Eq.(3.6) leads to the solution 



P 

V2 









\jk-^{y-yo) 



(3.9) 



In this case, we loose the soliton character of the profile. (The onset of the soliton 
depends on the parameters a and k). 
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Another interesting solution linked to the symmetry operator Vq arises for c 

3 
4' 



A; - 3/4 - \a\/2 > and a < 0, /of. Indeed, from Eqs. (3.6), (3.4) and (3.2) we 



have 

u 



V = -Lgf («-^)sn [y^{t -x),h], (3.10) 



where sn(-) is the sinus elliptic function of modulus h — 
Case b) By using the basis of invariants 

x = x, U{x) = ue-''^l^^ = «e-'*'/^^ 

W(x) = ve~''^/^^ = T;e-^*'/^^ (3.11) 

associated with the vector field V4, from Eqs. (1.1) we find the pair of reduced equa- 
tions 

i (u' + +kW+ (Qi|[/|2 + p\W\'^) U = 0, (3.12a) 

i (w' + +kU+ (a\W\^ + P\U\^) W = 0, (3.126) 

where U' ^ W ^ For a = /3 and assuming that k is real, Eqs. (3. 12) can 

be solved in the following way. First, let us divide (3.12a) and (3.12b) by U and W, 
respectively. Second, let us subtract the resulting equation corresponding to (3.12a) 
from that corresponding to (3.12b). Then, we obtain 

i {WW - UW) ^kiw"^ - U^) = 0. (3.13) 

By introducing now 

U=\{A + B), W=\{A-B) (3.14) 
into Eq.(3.13), we can determine B in terms of A, namely 

B = ^e-2^('=^+'5o), (3.15) 

where 5q is a constant of integration. Hence, putting the quantities (3.14) into 
Eq.(3.12a) and taking account of (3.15), after some manipulations we arrive at the 
solutions 

u = -^e*(2^-5o) ^iaa' COS ikx + 5o), (3.16a) 
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where a is a real constant. The physical role of these configurations remains to be 
explored. Here we limit ourselves to observe that the 'mass density' + |f p is of the 
Coulomb-type in the variable and is time independent. Moreover, the quantities 
(3.16) resemble the singular solution to the nonlinear Schrodinger equation expressed 
by formula (3.10) of ref.9. 

An interesting class of exact solutions related to the generator V4 can be found 
for q; 7^ /3 and k — iko, where ko is a real number. 

In this case, let us look for solutions to Eqs.(3.12) of the type 

U^pe'\ W^pe'^, (3.17) 

where p, 9 and 7 are real functions of x. Inserting (3.17) into Eqs.(3.12) yields 

o' 1 

!- + — + ko cosh -6) = 0, (3.18a) 
p 2x 

e' + kosm{-f-e) - (a + f3)p^^0, (3.186). 

7' - ko sin (7 - ^) - {a + p)p^ = 0, (3.18c) 

By integrating this system and using (3.11), we give the following pair of exact 
solutions to Eqs.(l.l): 

u^pe'^^ + ^\ V ^ ue'^''^^^^^^^''°^''°-''^\ (3.19) 

where 

/cosh2A;o(a;o — x) 



P 



X 



e = \ arcsin [cosh2A;o(xo - x)] + (« + /?) ^"^^^^0(^0 ^ (3.20) 

2 Jxi X 

and c, Xq, Xi, are real constants. 

At this point, let us deal with the special choice 9 = j. Then, Eqs.(3.18) lead to 
the solution 

e iL. i(a+/3)c^/ 7 dx' /n m\ 

u ^ c — ^e4xe ■'^1 ^ , (3-21) 

^/x 
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where c and xi are constants. (For xi — > oo, the integral at the r.h.s. of (3.21) 
becomes — Ei(2A;a;), where Ei(-) is the exponential integral function [ITUi). The 'mass 
density' corresponding to the solutions (3.20) and (3.21) is 



\u\' + \v\' = 2c''-^^^^^M^^ and \u\^ = \v\^ = c^'-^ . (3.22) 

X 



We point out that for A; = 0, the mass densities (3.22) behave as 1/x. Consequently, 
the presence of the parameter k = iko induces a change of the Coulomb-like mass 
density. 

Case c) A basis of invariants related to the symmetry operator Vi + V4 is given 

by 

7] = t - -X =t - -X , 

W{r]) = ^e-5*(''+^*') = t;e-^^(''+^^'). (3.23) 
By assuming a = P, and setting (3.23) into Eqs.(l.l), we get the reduced system 

U" - ^r]U + kW + a(|f/p + \W\'^)U = 0, (3.24a) 

W" - ^r]W + kU + a{\W\^ + \U\^)W = 0, (3.246). 

where U' = By requiring that U and W are real functions and U = W, 
Eqs.(3.24) lead to a special case of the second Painleve equation 



dz 



zi) + 2V^^ (3.25) 
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where z = 2i{^r] — k) and ip = 2^ ^J—a U. 

Case d) A set of invariants arising from the generator V5 is 

XX 

/= \/i{usm.kx + iv coskx) = \^{u sin kx + iv cos kx) , 

J= VS:{u cos kx — iv sin kx) = y/x{u cos kx — iv sin kx) , (3.26) 
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from which 

Equations (3.26) imply 

u — —j= (/ sin kx-\- J cos kx) , (3.27a) 

V X 

1 

V— {—I coskx + Jsinkx) . (3.276) 
Then, substitution from (3.27) into Eqs.(l.l) (for a — P) yields 

^(I- 2il') + ee'/' + AiH" + a + I J|2) / = 0, (3.28a) 

-Jt (J - 2^/) + e^V + 4eV" + a f + I J|2) J = 0, (3.286) 
2i ^ ^ 

where / = -^(C)> = ^' = % ^ ^ simple solution to the system 

(3.28) can be found supposing that /(^) and J{^) are real functions. In such a case, 

Eqs.(3.28) can be easily integrated. We have 



u 



-2 sin(A;x + 0) 



a t 



. -2 cos(/ca; + 0) 
v--^^- , (3.29) 

from (3.27), where is a constant and a may take both positive and negative values. 

4 Conclusions 

We have found some exact solutions to the system of equations (1.1) describing 
the propagation of waves in birefringcnt optical fibers. These equations have been 
analyzed in the framework of the Lie group theory. We have determined the associated 
symmetry algebra and the corresponding group transformations that leave Eqs.(l.l) 
invariant. Explicit configurations have been obtained both in the integrable and 
in the nonintegrable case. The subgroup (of the symmetry group) responsible for 
the sohton profile has been provided. Not all the configurations arising from the 
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symmetry approach have an evident physical meaning. This problem remains open. 
Notwithstanding, the knowledge of exact solutions to Eqs.(l.l) may be used with 
benefit as a guide for the development of perturbative techniques or for numerical 
calculations. Furthermore, the existence of exact solutions could be a challenge for 
trying to create new experimental patterns. 
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